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The atomic-to-molecular gas conversion is a critical step in the baryon cycle of galaxies, which sets the initial conditions for 
subsequent star formation and influences the multi-phase interstellar medium. We compiled a sample of 94 nearby galaxies with 
observations of multi-phase gas contents by utilizing public H1, CO, and optical IFU data from the MaNGA survey together 
with new FAST Hr observations. In agreement with previous results, our sample shows that the global molecular-to-atomic gas 
ratio (Rmo = log My,/My,) is correlated with the global stellar mass surface density u, with a Kendall’s v coefficient of 0.25 
and p < 107, less tightly but still correlated with stellar mass and NUV- r color, and not related to the specific star formation 
rate (SSFR). The cold gas distribution and kinematics inferred from the H1 and CO global profile asymmetry and shape do not 
significantly rely on Ry. Thanks to the availability of kpc-scale observations of MaNGA, we decompose galaxies into Hn, 
composite, and AGN-dominated regions by using the BPT diagrams. With increasing Ryo), the fraction of Hi regions within 
1.5 effective radius decreases slightly; the density distribution in the spatially resolved BPT diagram also changes significantly, 
suggesting changes in metallicity and ionization states. Galaxies with high Rmo tend to have high oxygen abundance, both at one 
effective radius with a Kendall’s r coefficient of 0.37 (p < 107°) and their central regions. Among all parameters investigated 
here, the oxygen abundance at one effective radius has the strongest relation with global Ri. The dependence of gas conversion 
on gas distribution and galaxy ionization states is weak. In contrast, the observed positive relation between oxygen abundance 
(u) and Roi indicates that the gas conversion is efficient in regions of high metallicity (density). 
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1 Introduction 


Gas is the raw material that forms stars in galaxies, which 
plays a crucial role in galaxy formation and evolution (e.g., 
[56, 57, 66,69, 94]). The cold gas predominantly consists of 
neutral atomic (H1) and molecular hydrogen (H2). Within 
the baryon cycle, atomic gas collapses to form dense molec- 
ular gas, which directly fuels star formation in galaxies. 
The tight relations between star formation and molecular 
gas, in terms of surface density (Schmidt-Kennicutt law, 
[38, 56, 66, 70, 83, 104 ]) and global mass (e.g, [43, 97]), have 
been extensively studied, but insufficient research has been 
conducted to adequately investigate the connection between 
the atomic-to-molecular gas conversion process (H r-H5) and 
factors such as metallicity or ionization states. 

The H i-H? conversion primarily depends on the mid-plane 
pressure [11,12,41], radiation field [41], and metallicity [62]. 
The dominant factor is the mid-plane pressure, which is pro- 
portional to the surface density (both gas and star) after as- 
suming a thin disk with balanced gas and stars [11]. The 
radiation field influences the molecular gas dissociation, and 
metals/dust act as the catalyst in the H 1-H; conversion. Leroy 
et al. [66] found that the surface density ratio of H5 and H1 
for spiral and dwarf galaxies increases with the increasing 
of local stellar mass surface density and mid-plane pressure, 
while it decreases with the increasing of galactic radius. If the 
Hi surface density exceeds ~ 9 Mo pc? at solar metallicity, 
molecular hydrogen can stably exist against self-shielding, as 
supported by both theoretical [62, 63] and observational [9] 
evidence. Moreover, both semi-analytic models [62, 87] and 
hydro-dynamical simulations [64] show consistent results. 

The global values of molecular-to-atomic gas ratio (Rio = 
log My,/My,) describe the efficiency of gas conversion and 
constrain the physical properties of the interstellar medium 
(ISM, [15p. Aga shows a weakly increasing trend with 
the global stellar mass M., stellar mass surface density us, 
and NUV-r color for galaxies with log (M./Mo) = 10.0 [95]. 
These relations become stronger while extending to the low 
mass end and with a larger sample [15, 16, 23]. Stark et al. 
[110] found a positive relation between Rmo; and enhance- 
ments of central star formation in spirals, but the enhance- 
ment is not applicable for massive early-type galaxies. If 
Ls = 1097 Mo pc^?, the cold gas reservoir is depleted or dis- 
pelled; otherwise, Rmo is positively correlated with u, [94]. 
Additionally, the spatial distribution of H 1 and H5 may play 
a crucial role in their conversion. 

The optical diagnostics inferring galaxy ionization states 
and metallicity may play an important role in understanding 
the physics of gas conversion. Ionization states of gas in a 
galaxy can be inferred through optical diagnostics such as 
the BPT diagrams [6]. These diagnostic tools characterize 
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the ionization states by using strong optical emission lines, 
such as H œ A 6564, H B A 4862, [S n] 44 6716, 6731, 
[N u] AA 6548, 6583, and [O m] 44 4959, 5007. For in- 
stance, the line emission from galaxies or regions can be 
categorized into Hu, composite, and AGN-dominant emis- 
sions based on their positions in the [O m]/H versus [N 
iHa ([N u] BPT) plane following Kewley et al. [58] and 
Kauffmann et al. [55] or using the P1P2 BPT diagram fol- 
lowing Ji & Yan [52]. Moreover, the [S m]-based BPT ([O 
m]/H8 versus [S n]/Ho) classifies galaxies or regions into 
Hu, Low-Ionization Nuclear Emission Line Regions (LIN- 
ERs), and Seyfert [59]. In some special cases, shocks can 
show up even in the star-forming regions of the BPT diagrams 
[2]. Typically, H u emission-dominated galaxies or regions 
primarily exhibit emission from massive young stars, while 
AGN-dominant regions are driven by AGN photoionization 
or shock, and composite components are results of a mixture 
of star formation, shock excitation, and/or AGN activity [60]. 

Theoretical models take metallicity as a controlling param- 
eter in H» formation [41, 62, 87], which is supported by ob- 
servations in nearby galaxies. Rmoy and gas metallicity are 
tightly correlated [15], because metallicity serves as an indi- 
cator of dust content, which facilitates the H 1-H» conversion 
as a catalyst [62]. At a given stellar mass, a high star forma- 
tion rate (SFR) is associated with low metallicity [28, 39] but 
a high both atomic and molecular gas fraction [96]. Thus the 
complex effects of gas content and metallicity on the atomic- 
to-molecular gas conversion are non-trivial. The influence of 
the metallicity distribution on the H i-H» transition remains 
unknown. So more studies are critical to reveal the connec- 
tion between gas conversion and metallicity. 

Detailed studies on this field require high-resolution op- 
tical integral field unit (IFU) data and cold gas data with 
a large sample size. The Mapping Nearby Galaxies at 
APO (MaNGA, [1, 19]) survey and nearby cold gas sur- 
veys provide spatially resolved metallicity and gas content 
information for nearby galaxies. The MaNGA project uti- 
lizes IFU spectroscopy to map the detailed composition 
and kinematic structures of 10,010 nearby galaxies. Its 
H1 follow-up (H1-MaNGA: [76, 111]) collects 3669 unique 
Hı observations from GBT and Arecibo. Moreover, the 
ALMA MaNGA Quenching & Star-Formation Survey (AL- 
MaQUuEST) utilizes high-resolution spatially resolved opti- 
cal spectroscopy from MaNGA and CO (J = 1-0) follow-up 
observations with ALMA in the C43-2 configuration to study 
stellar and gas properties on the same physical scales [70,71]. 
The IFU and CO observations of MaNGA and ALMA have 
a beam size of 2.5" (1” ~ 1 kpc) with 0.02x z <0.05 and a 
field of view of 50". For galaxies with MaNGA IFU, CO, and 
H1 data, we could investigate the relation between gas con- 
version and optical diagnostics. However, the ALMaQUEST 
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survey lacks complete and homogeneous H1 observations. 

In this paper, we present new observations of Hı global 
spectra for the 37 galaxies in the ALMaQUEST survey using 
the Five-hundred-meter Aperture Spherical Telescope (FAST, 
[81], [68]). To build a large sample for statistical analy- 
sis, we match the MaNGA survey with ALMaQUEST, the 
CO follow-up of the GALEX Arecibo SDSS survey (xCOLD 
GASS: [95, 97]), the APEX Low-redshift Legacy Survey 
for MOlecular Gas (ALLSMOG: [16, 25]), and “JCMT dust 
and gas In Nearby Galaxies Legacy Exploration” (JINGLE: 
[32,98]) surveys to investigate the dependence of Rmoi on gas 
distribution and metallicity. In Section 2, we introduce the 
sample properties and the cross-match results. The Hı data 
processing and spectral measurements are presented in Sec- 
tion 3. We investigate the dependence of molecular-to-atomic 
gas ratio in Section 4. The main results are summarized in 
Section 5. This work adopts the following parameter for a 
ACDM cosmology: Q,,- 0.315, Qa= 0.685, and Ho = 67.4 
km s^! Mpc"! [86]. We adopt a Chabrier initial mass function 
(IMF, [24]) to uniformly derive SFRs and stellar masses. 


2 TheSample 


We compile a sample of galaxies with atomic and molec- 
ular gas measurements from Hı and CO surveys to con- 
struct a comprehensive sample for investigating the atomic- 
to-molecular gas conversion (H r-H5). The molecular gas in 
nearby galaxies has been extensively studied through various 
large observation programs such as the ALMaQUEST sur- 
vey [70,71], xCOLD GASS [95,97], ALLSMOG [16, 25], 
and the JINGLE survey [32,98]. The xCOLD GASS sam- 
ple covers a stellar mass range of 10°° — 10!!?^ Mọ with a 
redshift of 0.01—0.05 and includes star-forming, green valley, 
and quenched galaxies. The H 1 surveys primarily consist of 
xGASS [22,23] and H -MaNGA [76,111], which combines 
new GBT observations with existing data from the Arecibo 
Legacy Fast ALFA (ALFALFA) survey [48,49]. We collect 
the H 1 and CO data of our sample galaxies from the datasets 
listed above. 


2. The ALMaQUEST and Supplementary Samples 


We obtained the FAST Hı spectral line observations for 37 
galaxies, which are selected from 51 galaxies from the AL- 
MaQUEST survey [71] and supplemented with 2 galaxies 
(MaNGA 8439-6102 and MaNGA 8250-6104) from Gao et 
al. [44]. These 51 galaxies include 46 published systems 
in the original ALMaQUEST Survey and 5 additional galax- 
ies from the extended ALMaQUEST merger sample [115], 
the latter of which are either not covered by previous H 1 ob- 
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servations or have a detection of signal-to-noise ratio (S/N) 
less than 3. The two galaxies included in Gao et al. [44] 
exhibit a gas-rich nature in their molecular phase (CO 2- 
1) and are actively star-forming. Their stellar masses and 
redshift fall within similar ranges (10.0< log M,/Mo«11.5, 
z ~0.04) compared to the ALMaQUEST sample. We also 
include these two systems in our analysis. 

These 37 galaxies were not included by H -MaNGA data 
release 1 [76] but fall within the declination range of — 14? to 
66^. This selected sample (referred to as the ALMaQUEST 
sample in this paper) consists of galaxies with a stellar mass 
range of 10.0 x log (M./Mo) € 11.5 and a redshift range of 
0.02 x z € 0.13. This sample primarily consists of galaxies in 
the green valley, star-forming main sequence, and star-burst 
categories [40, 70-72]. We adopt the values of stellar masses 
and SFRs by using the Chabrier IMF instead of Salpeter IMF 
used in previous work through utilizing the conversion of 
Madau & Dichinson [73]. 

To investigate the dependence of H 1-H) conversion, we 
also compile a set of galaxies with similar properties to the 
ALMaQUEST sample. In addition to the ALMaQUEST sur- 
vey, galaxies with both H1 and CO observation are mainly 
from the xCOLD GASS, ALLSMOG, and JINGLE surveys. 
Sample galaxies selected from these surveys are with rela- 
tively uniform sample selection criteria, accordant global H1 
and CO profiles, and a large sample size that covers a stellar 
mass range of 9.0 x log (M./Mo) x 11.5. Among these sur- 
veys, the xCOLD GASS survey provides both accessible H1 
and CO spectra. 

We adopt a constant CO(1-0) to Hz conversion factor: 
aco = 4.35 Mo (Kkms^! por following Bolatto etal. [14] 
and references therein. Our sample has an oxygen abundance 
12+log (O/H) 2 8.4 (see details in Section 4.4), thus the value 
of aco does not depend on galaxy metallicity [14, 103]. For 
galaxies in the ALLSMOG and JINGLE surveys that only 
have CO(2-1) observations, we convert the intrinsic bright- 
ness luminosity Lcoo.1) to Lco(-0) by assuming a ratio of r21 
= Lcoo-)/Lcou-oj7 0.8 as suggested by Cicone et al. [25] 
for nearby normal star-forming galaxies, and the recent spa- 
tially resolved studies from Leroy et al. [67] also give similar 
but slightly lower values. We investigate the dependence of 
cold gas conversion on galaxy ionization states (Section 4.3) 
and metallicity (Section 4.4). 
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2.2 Sample Definitions 


In the following discussions, we refer to galaxies with both 
published Hr and CO spectra as the “HICO-Spec” sample. 
This sample comprises galaxies with H1 and CO observa- 
tions (with at least one detection in Hi or CO) from the 
ALMaQUEST and xCOLD GASS samples. After remov- 
ing 2 duplicate galaxies, the HICO-Spec sample contains 412 
galaxies with both H1 and CO observation, of which 310 
galaxies have both Hı and CO detections. In Section 4.2, 
we investigate the dependence of H r-H» conversion on inte- 
grated profile shape and asymmetry using this sample. 


xGASS 
*  ALMaQUEST detection 
2 yz  ALMaQUEST non-detection 
+ HICO-Spec — 
1 +  HICO-MaNGA— — — —e—-$— $- = x. 
E UO D- 
" -- 
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© -r AD = 
= —1 = mx E 
ac - >. 
LL. t F : 
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9 x » 
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10. 110 11.5 
log (M «/Mg) 
Figure 1 The distribution of galaxies in the HICO-MaNGA (red) and 


HICO-Spec (black) samples in the space of log SFR vs. log M.. The de- 
tection and non-detection of the ALMaQUEST galaxies are shown as filled 
green and open cyan stars, respectively. In comparison, we include MaNGA 
(blue scale and contours) and xGASS (small grey dots) galaxies in the plot. 
The black solid and dashed curves show the star-forming main sequence [94] 
and its deviation at + 0.4 dex. 


On the other hand, we refer to the “HICO-MaNGA” sam- 
ple as the cross-match result of the cold gas surveys (AL- 
MaQUEST, xCOLD GASS, ALLSMOG, and JINGLE) and 
the MaNGA survey [99,100]. During the cross-matching pro- 
cess, we utilize TOPCAT [114] and require a maximum pro- 
jected optical position separation of 3" and a velocity sep- 
aration of 500 km s~!. We compile the sources in the cold 
gas catalogs and remove 4 duplicated galaxies, resulting in 
a total of 390 galaxies with both H 1 and CO detections. We 
cross-match these 390 galaxies with the MaNGA catalog of 
10,010 unique galaxies by their positions and redshifts. This 
step yields 32 galaxies from xCOLD GASS, 7 galaxies from 
ALLSMOG, 38 galaxies from ALMaQUEST, and 17 galax- 
ies from the JINGLE survey. Thus the HICO-MaNGA sam- 


January (2024) Vol. xxx No. xx 


000000-4 


ple consists of 94 galaxies with both H1 and CO detection 
as well as optical IFU data from MaNGA. Figure 1 shows 
the sample distribution in the space of log SFR vs. log M. 
for the HICO-Spec and HICO-MaNGA samples. The sample 
size of HICO-Spec is slightly larger than that in Saintonge et 
al. [95]. The galaxy number of the HICO-MaNGA sample 
is 10% of the xCOLD GASS survey [23], but all galaxies in 
the HICO-MaNGA sample have MaNGA IFU observations. 
Meanwhile, our sample covers a similar dynamical range in 
stellar mass and SFR compared to the star-forming galaxies 
in the MaNGA survey. 


2.3 FAST H Observations 


In addition to the existing Hı observations from GBT and 
Arecibo, we obtained FAST Hı spectral line observations 
for the ALMaQUEST and Supplementary Samples. The 
observation is designed to reach a detection threshold of 
log Mu,/M. z 1% following the strategy of GASS and 
xGASS [22, 23]. In 2020, we conducted the initial Hı 
follow-up observation using FAST (PI: Zheng, project code: 
PT2020.0102). We performed additional observations on 
two galaxies (MaNGA 7977-12705, MaNGA 8077-6104) 
that were previously detected by H -MaNGA (GBT) to ver- 
ify the calibration procedure. It should be noted that a portion 
of the data was affected by strong radio frequency interfer- 
ence (RFI), primarily caused by the refrigerating dewar in the 
compressor. However, this issue has been partially resolved 
in 2021 [128]. Taking into account the data obtained in 2020, 
we have adjusted our observation strategy and proposed a 
new FAST observation in 2022 led by PI Zheng (project code: 
PT2022. 0091). 

The observations were conducted using the entire wide 
band of 500 MHz, employing the on-off mode and the 19 
beam receiver at the L band (1.05-1.45 GHz). Each beam 
has a size of 2.9' and a sampling time of approximately 0.1 
seconds. During each second, the 10 K (high CAL) noise 
diode is on for the first «0.1 seconds and off in the remain- 
ing 0.9 seconds. The spectrum encompasses 65,536 channels 
and has a spectral resolution of 7.6 kHz (~1.6 km s^!) for 
the dual-polarization observations. The separation between 
the center of beam M01 and one of the outer beams M14 is 
11.6’. 

In each cycle of observation, the central beam MOI points 
on the source and the outer beam M14 points away from the 
source in the “source on" mode. Conversely, the beam M14 
points on the source and the beam MOI points away from 
the source in the "source off" mode. The switching time be- 
tween source on and source off modes is 30 seconds. Be- 
sides the switching or overhead time, the observation beam 
keeps pointing on the source, effectively doubling the inte- 
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Table 1 Log of FAST Hı Observations for 37 Galaxies 


Plate-IFU R.A. Decl. z Dij  HrMaNGA Sample Tin Hı Flag 
(J2000) (J2000) (Mpc) (seconds) 

(1) (2) (3) (4) (5) (6) (7) (8) (9) 
7815-12705 = 21:15:57.71 = +09:32:35.1 — 0.02955 134 0 0 120x3+120x2 1 
7977-3704 22:11:11.70  «11:48:02.6 0.02724 123 2 0 120x5+120x3 1 
7977-9101 22:04:29.50  +12:26:33.5 0.02656 120 2 0 120x1+120x3 1 
7977-12705 22:11:34.28 «11:47:45.3 0.02724 123 1 0 120x1+120x1 1 
8077-6104 02:48:07.87  -00:45:08.3 0.04601 211 1 0 180x9+180x4 1 
8077-9101 02:46:34.35 | —00:50:36.7 0.04323 198 2 0 180x8+120x4 1 
8081-3704 03:19:17.15 — —00:58:10.7 — 0.05400 249 0 0 180x8+180x4 1 
8081-6102 03:19:45.63 — —00:04:37.9 — 0.03719 169 2 0 120x3 3 
8081-9101 03:11:05.32 —00:32:47.5 0.02846 129 2 0 120x4+120x4 1 
8081-12703 . 03:21:33.93  -00:10:42.1 0.02558 116 2 0 120x7+120x4 1 
8082-12704  03:19:47.89  —00:13:16.1 0.13214 642 0 0 180x8 3 
8083-12702 . 03:20:58.90 | —00:22:03.7 0.02104 95 0 0 120x1+120x2 1 
8084-3702 03:22:32.79  —00:00:04.4 0.02206 99 2 0 180x9+180x3 1 
8084-6103 03:22:58.00  +00:03:14.9 0.03593 164 2 0 180x8+180x10 2 
8155-6101 03:35:15.39  -01:13:43.0 0.03740 170 1 0 120x2+120x3 1 
8155-6102 03:30:29.13 +00:45:07.4 0.03081 140 0 0 120x5 1 
8156-3701 03:42:22.15  —00:34:59.5 0.05273 243 0 0 180x9+180x10 1 
8241-3703 08:25:50.69  +18:10:00.1 0.02911 132 0 0 180x10 1 
8241-3704 08:26:16.54  +17:21:44.8 0.06617 308 0 0 180x8+180x4 1 
8615-3703 21:23:18.34  +401:15:17.9 0.01845 83 0 0 180x8 1 
8615-12702  21:20:38.27 = +01:02:50.2 0.02095 94 0 0 120x5+120x3 1 
8616-6104 21:31:55.33 400:12:49.6 0.05426 250 0 0 180x9 1 
8618-9102 21:17:05.15 — +09:58:20.3 0.04334 198 1 0 180x10+180x10 1 
8623-6104 20:47:07.44 — 400:18:01.7 0.09704 461 0 0 180x9 3 
8623-12702 . 20:40:52.10 . 400:39:10.1 0.02691 122 0 0 120x3 1 
8655-3701 23:47:00.44  —00:26:50.6 — 0.07149 334 0 0 180x9--180x3 1 
8655-9102 23:52:53.26  —00:22:56.8 0.04505 206 2 0 180x10+180x10 1 
8728-3701 03:50:47.77 -07:01:43.6 0.02833 128 2 0 120x2+120x4 1 
8950-12705 12:58:55.95 +27:50:00.4 0.02528 114 2 0 120x2 3 
8952-12701 . 13:38:44.12 +26:19:42.7 0.02856 129 2 0 120x2 3 
9194-3702 03:08:07.07  +00:27:22.4 0.07454 349 0 1 180x10+180x10 1 
9195-3702 01:51:22.27  413:03:37.2 0.06434 299 0 1 120x4+180x4 2 
9512-3704 09:12:36.30  +00:18:34.2 0.05463 252 2 1 180x9 3 
8153-12702  02:41:12.96  -00:52:37.3 0.03823 174 1 1 180x8+180x10 1 
7975-6104 21:39:33.98 +10:29:00.5 0.07920 372 0 1 180x10 3 
8439-6102 09:31:06.76  +449:04:47.1 0.03415 155 2 2 120x2+180x10 1 
8250-6104 09:21:38.74  +443:43:34.1 0.04008 183 2 2 120x7+180x4 1 


Column (1): the MaNGA Plate-IFU. Columns (2)-(3): equatorial coordinates (J2000). Column (4): redshift from SDSS. 
Column (5): luminosity distance by adopting cosmology of the Planck Collaboration [86]. Column (6): the status of H1 obser- 
vation in H1-MaNGA - 0: not observed, 1: detection, 2: non-detection. Column (7): the original sample — 0: ALMaQUEST, 
1: the additional sample in ALMaQUEST, 2: Gao et al. [44]. Column (8): the total integration time in unit of second of each 
cycle with the numbers showing the integration time per cycle and total cycles. The initial and subsequent ones indicate the 
observation times of 2020 and 2022, respectively. Column (9): the combined quality of FAST H1 observation — 1: detection; 


2: non-detection, 3: contaminated by RFI. 
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gration time compared to observations without pointing on 
and off simultaneously. Previous studies [134, 137] have also 
employed the use of other beams as off-target points to save 
observation time. 

We process the FAST data to obtain the final global spec- 
tra. To mitigate the effects of RFI, we remove values with a 
deviation from the median value larger than 1.50 per chan- 
nel in each observation cycle. This pre-processing step was 
necessary for the FAST data of 2020. We then interpolate 
the remaining FAST raw data using the values of the nearest 
neighbors to supplement data points. We calibrate the flux 
intensity using data from the on-off noise diode and apply- 
ing aperture efficiency correction taking the beam, frequency, 
and position into consideration [53,54]. Then the data of the 
first two polarizations is averaged to derive the preliminary 
spectrum. 

To remove the standing wave), we fit a sine function after 
masking the RFI and the velocity range of the potential emis- 
sion line signal. If the removal of the standing wave did not 
decrease the noise level of a specific spectrum, we did not 
adopt the standing wave subtraction. Subsequently, a poly- 
nomial of order 1-3 was used to fit the spectrum and flatten 
the baseline. The degree of polynomial order is determined 
by selecting the degree with the minimum Bayesian Infor- 
mation Criterion (BIC, [107]), which balances the goodness 
of the fit (reduced chi-square) with the number of parame- 
ters in the polynomial fitting. We convert the frequency to 
heliocentric velocity by using the optical velocity convention 
and performing the Doppler correction. For each source, we 
stack the spectrum of each cycle (at least two cycles for each 
galaxy by using M01 on and M14 on) by weighting the spec- 
trum according to its noise level [18,42]. 

We observed 37 galaxies with a total observation time of 
45 hours, which includes the source on, source off, and over- 
head time. The integration time for each galaxy ranges from 2 
minutes to 1 hour. Table 1 provides the basic information and 
observation details for these galaxies. For a RFI masked and 
baseline subtracted profile, we search for the emission line 
signal by rebinning the spectra to a channel width of ~ 20 km 
s^!. The noise level of the spectrum (Ospec) is the standard 
deviation obtained from the best-fit half-Gaussian fitting to 
the negative values of all flux intensities within a range of + 
500 km s^! around the optical central velocity. If the rebinned 
spectrum contains at least three channels with flux intensities 
greater than 30 gpec, it is considered a detection; otherwise, it 
is classified as a non-detection. The FAST Hı spectra quality 
is listed in Column (9) of Table 1. 

The new FAST observation yields 27 detections, 1 absorp- 
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tion signal, and 2 non-detections, and all thses spectra are 
public ?. Among the 53 galaxies in the ALMaQUEST sam- 
ple, we have acquired 42 emission line signals, 1 absorption 
signal, 7 non-detections, and 3 RFI-contaminated ones by 
combining the data of H -MaNGA [76,111] and new FAST 
observations. There is no useful signal for MaNGA 8082- 
12704, MaNGA 8623-6104, and MaNGA 7975-6104 due to 
heavily RFI contamination, and thus were excluded from fur- 
ther discussion. The distribution of ALMaQUEST detection 
and non-detection of H i signal is shown in Figure 1. Around 
half of the non-detections are massive green-valley galaxies 
with log M,> 10.8, which may suggest that galaxy quenching 
is due to gas depletion. 


3 Results and Analysis 


3.1 FAST Hı Spectra 


The final results of FAST observations are shown as blue pro- 
files in Figure 2 and 3, which includes emission, absorption, 
and non-detection. To align with the spectra obtained from 
Arecibo and GBT, we rebin the spectra to a channel width 
of ~ 6 km s !. The measurements of spectral noise level o 
(in units of mJy per channel) follows the same approach as 
Ospec, but with a narrower channel width AV ~ 6.4 km gt 
The black profiles in Figure 2 are observations from the H 1- 
MaNGA survey. The comparison of the blue and black global 
profiles in Figure 2 demonstrates the excellent consistency 
between the FAST and archival (Arecibo and GBT) spectra, 
but the new data exhibit a relatively lower noise level. 


While the H -MaNGA data indicate a non-detection for 
MaNGA 8155-6101, the FAST data reveals a high S/N Hı 
absorption signal by reducing the noise level by half, which is 
submerged in the noise in the GBT spectrum (we will return 
to a discussion of this absorption signal in Section 3.3). The 
discrepancy in spectra between FAST and GBT demonstrates 
the capability of FAST in detecting new signals. Regarding 
MaNGA 8084-6103, the new FAST observation further con- 
strains the upper limit of the H 1 mass to My rig = 1085 Mo, 
which is 1.1 dex deeper than that of H - MaNGA. The upper 
limit of H1 mass is calculated by assuming a line width of 200 
km s^! and a 3c detection following H -MaNGA [76,111]. 
So the upper limits for the H 1 mass is 


D 
Muy 11im/Mo = 2.356 x 10° x (— yx 
Mpc 


1) a kind of periodic sinusoidal feature appeared in the baselines, which is caused by reflected broadband signals in the focal structure of the telescope 


entering the receiver system with time delay [17]. 
2) https://github.com/NiankunYu/ALMaQUEST.FAST-HI2024 
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Figure 2 The Hı spectra of galaxies with both FAST (blue) and H1-MaNGA (black) observations. In each panel, the MaNGA IFU plate ID is shown in the 
upper left corner, the optical central velocity is shown as the red vertical line. Among these profiles, the spectrum of MaNGA 8084-6103 shows no detection 
above 3-o limit, MaNGA 8155-6101 shows an absorption line signal, and the Ht emission of MaNGA 8728-3701 is contaminated by its companion. The 
optical centers of the companion galaxies of MaNGA 8155-6101 and MaNGA 8728-3701 are shown as green vertical lines. 
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identified as a non-detection. 


200 x (AV/km s~’) x 3or/(1000mJy) 
(Cim 
y kms 


and the H 1 mass upper limit of 7 non-detections are shown 
in Table 3. Additionally, the Hı profile of MaNGA 8728- 
3701 from H -MaNGA is classified as a non-detection, while 
the strong H r emission at 8400-8600 km s^! originates from 
its companion, MaNGA 8728-12701. The new FAST data 
provides further confirmation that the weak emission around 
8300 km s^! may be from MaNGA 8728-3701 and give a 


conservative H 1 mass upper limit of 10?? Mo. Due to the 


inclusion of the emission at 8200-8600 km s^! , our H 1 mass 
upper limit exceeds that of H -MaNGA. 


3.2 Measurements of Global Profiles 


The H 1 emission line spectra of the new FAST supplements 
and xCOLD GASS datasets are uniformly measured using 
the curve-of-growth method [131, 133]. This method defines 
the profile center as the flux-intensity weighted velocity, and 
constructs the curve of growth by accumulating flux inten- 
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Table 2 Parameters Derived from the H1 Spectra of the ALMaQUEST Sample 

Plate-IFU Ve F Vas Ar Ac Cy K S/N c log Mui 

(kms!) (ykms7!) (kms?) (mJy) (Mo) 
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) 
7815-12705 8785 x3  2.10:0.33 329 x20 1.094+0.08 1.30+0.13 2.95€0.21  -0.04x0.02 108.9 0.2 9.9 € 0.1 
7971-3704 8114x6 0.50+0.08 285417 1.04+0.08 1.03+0.12 2.932023  -0.04 + 0.02 182 0.4 9.2 x 0.1 
7977-9101 788849 0.294005 348426 1.00+0.09 1454017 2.4940.22  -0.07 x0.02 13.1 0.3 9.0+0.1 
7977-12705 8094+3 46140.73 371422 1.00+0.07 1.0040.10 3.06+0.22 -0.03 + 0.02 66.1 09 10.2 0.1 
8077-6104 13724+3 1.27 0.20 153+9 1.05+0.08 1.08+0.11 2.92+0.21 -0.03 + 0.02 73.3 0.2 10.1 0.1 
8077-9101 12849+7 0.50+0.08 438+33 1.004+0.08 1.25+0.15 2.62+0.25  -0.05 + 0.03 20.8 0.3 9.6 + 0.1 
8078-6103 8516+10 2.0340.33 276 x24 1.144013 1.044013 422+40.39 0.04 + 0.02 10.4 3:3 9.9 + 0.1 
8078-12701 7995+11 341+053 408425 1.08+0.10 1.36x023  Á 2.332020 -0.10+0.02 10.6 4. 10.0+0.1 
8081-3704 16160+3  Á0.75x0.12 190412 1.19+0.09 1.324015 4.83 +0.34 0.09 + 0.02 43.2 03 10.0+0.1 
8081-9101 845246 0.914014 393425 1.124008 1.104012 2.814021 -0.04+0.02 24.2 0.5 9.5 x 0.1 
8081-9102 1013148 240+0.38 390424 1.02+0.08 1.12+0.13 2.69x021  -0.05x0.02 18.5 18 10.1401 
8081-12703 7578+6 0.334005 310419 1.00+0.08 1.024014 2.4940.20  -0.08 + 0.02 17.4 0.3 9.0+0.1 
8082-6103 720745 142+0.22 268419 1.3040.11 1.744020 4.2440.39 0.07 + 0.02 14.5 1.8 9.6+0.1 
8082-12701 803449 2.884045 302421 1.104009 1.174015 2.80+0.25 -0.05 + 0.02 11.5 41 10.0+0.1 
8083-6101 793846 3.724059 382426 1.0140.08 1.384015 2.9640.23 -0.03 + 0.02 24.6 23 10.140.1 
8083-9101 11445413 1.54€ 024 357430 1.1640.15 1.004015  2.99x0.39 -0.05 + 0.03 11.9 18 10.0+0.1 
8083-12702 6312+3 43640.69 183411 1.064007 1.174012 2.954021 -0.0340.02 168.2 04 10.0+0.1 
8084-3702 6562+4 1.184019 261416 1.05+0.08 1.0640.11 3.0840.22 -0.02 + 0.02 42.3 0.4 9.4 x 0.1 
8084-12705 754346 2.884045 405426 1.18+0.09 1.324014 3.00023 -0.01 + 0.02 24.4 1.7 9.9+0.1 
8086-9101 11979 x13 1.184018 365431 À 1.00-0.00 1.3540.25 4.03 + 0.36 0.02 + 0.02 12.5 1.2 9.9+0.1 
8155-6102 9165+4 2.054032 167410 1.004+0.07 1.014011 3.0240.22  -0.03 + 0.02 33.7 10 10.0+0.1 
8156-3701 1571244 1.224019 259416 1.3440.10 1.724018 5.22+0.37 0.10 + 0.02 67.6 0.2 10.2+0.1 
8241-3703 8655+3 2.25 + 0.36 13548  À1.07x0.08 1.04+40.11 2.97+0.21 -0.03 + 0.02 54.7 0.7 10.0+0.1 
8241-3704 1974844 1362x022 293418 1.10+0.08 1.154012 3.66+0.26 0.01 + 0.02 56.6 03 10.5+0.1 
8450-6102 1253347 0.78+0.14 101413 1.0940.19 1.2940.33 5.27 = 0.90 0.08 + 0.03 8.6 2:3 9.8 x 0.1 
8615-3703 5483 +4 1.214019 248415 1.02+0.07 1.01+0.10 3.924028 0.03 + 0.02 50.8 0.3 9.3 € 0.1 
8615-9101 9961 £6 1.88+0.29 402-24 1.044+0.08 1.064012 2.932022  -0.05x0.02 20.0 13 10.0+0.1 
8615-12702 6198 £3 7.5541.19 257415  1.012z0.07 1.0240.10 3.0740.22 -0.0440.02 303.8 04  102z0.1 
8616-6104 1618927 0.51€ 0.08 408425 1.05+0.08 1.12x0.13 2.8440.21 —0.04 + 0.02 21.8 0.3 9.9 + 0.1 
8616-9102 904844 3.234051 212413  Á 1L012z0.07 1.004010 3.032022  -0.03x0.02 36.0 16 10.2+0.1 
8616-12702 913146 2.964047 450428 1.22+0.09 1.83+40.20  Á 2.842023 -0.02 + 0.02 26.5 16 10.1401 
8618-9102 1291325 040+0.06 273416 1.00+0.07 1.004010 3.44+40.25 -0.01 + 0.02 26.4 0.2 9.6 x 0.1 
8623-12702 8042+3 5.344084 391423  1.07x0.07  1.25x0.13  2.87x0.20 -0.0640.02 156.4 04  Á 103201 
8655-3701 2133147  0.27x0.00 216416  1.00-0.00 1.004012 4594041 0.08 + 0.02 11.1 0.4 9.8 x 0.1 
8655-9102 1343745 0.30+0.05 207413 1.06+0.08 1.154013 3.122024  -0.02 + 0.02 18.6 0.2 9.5 x 0.1 
8655-12705 1341329  Á1.0720.17 412428 1.29+0.11 1.304017 3.312029 0.05 + 0.02 13.3 12 10.0+0.1 
8952-6104 8460+4 1.70+0.27 62+5 1.184016 1.304019 3.67+0.46 0.03 + 0.02 14.3 2.9 9.8+0.1 
9194-3702 22245+4 0.704011 264417 1.1440.08 1.224013 3.92+0.30 0.05 + 0.02 48.2 0.2 1032401 
8153-12702 1143524 040+0.06 184411 1.01+0.07 1.044011 3.132023 -0.02 + 0.02 33.1 0.2 9.4 x 0.1 
8439-6102 10145x5 0.30+0.05 306418 1.05+0.08 1.394015 3.9840.31 -0.00 + 0.02 26.9 0.2 92 x 0.1 
8250-6104 11957 x7 0.59+0.09 376424 1.00+0.08 1.04+0.11 3.56 = 0.27 0.01 + 0.02 22.9 0.3 9.7 x 0.1 


Column (1): the MaNGA Plate-IFU. Column (2): Flux intensity-weighted central velocity. Column (3): Total global flux of the H1 line. Column (4): 
Velocity width measured at 85% of the total line flux. Column (5): Global flux asymmetry. Column (6): Flux distribution asymmetry. Column (7): 
Profile concentration. Column (8): Profile shape. Column (9): S/N of the profile. Column (10): Noise level of the profile at a channel width of ~6.4 
km s7!. Column (11): H1 mass and its uncertainty; for a distance uncertainty of 10% and a flux uncertainty of 15%, the typical uncertainty of log My, 


is 0.1 dex. 
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sities from the center to two sides, then quantifies the total 
flux, line width, asymmetry, and profile shape. We remove 
two galaxies: MaNGA 8728-3701 and MaNGA 8155-6101 
due to contamination of companions and absorption signals 
in section 4. 


The total flux F (in units of Jy km s^!) is determined as 
the median value of the flat part of the growth curve. The line 
width is defined as the velocity width that encloses a char- 
acteristic fraction of the total flux. For instance, line widths 
V25 and Vas correspond to the velocity widths enclosing 25% 
and 85% of the total flux, respectively. The profile asymme- 
try is quantified as the ratio of the integrated fluxes for two 
sides (Ar 21.0) and the ratio of the slopes of the rising part 
of the curve of growth for two sides (Ac 21.0). We quantify 
the emission line profile shape as the degree of the concen- 
tration of the line profile Cy and the integrated area between 
the normalized curve of growth and the diagonal line of unity 
K. With the increasing of Cy or K, the profile transits from a 
double-horned profile to a single-peaked profile. The S/N of 
a spectrum is calculated as 


1000F 
syn = 100 


-——. Q) 
c AV N2N 


We apply a second-order polynomial correction to account 
for systematic bias in the total flux, line width, profile asym- 
metry, and profile shape when the spectrum has a S/N below 
40 following Yu et al. [133]. The H 1 mass is calculated as 
My, = 2.356 x 10D? F Mo [91], where Dz is the luminos- 
ity distance with a unit of Mpc. The FAST beam roughly 
corresponds to a physical scale of 200 kpc, which is enough 
to cover all Hr emissions within the host galaxy, thus the 
aperture correction is not necessary for our targets. The fi- 
nal measured parameters of all detections obtained from our 
Hı spectra are presented in Table 2. 


The H1 mass upper limits of the ALMaQUEST sample 
cover a range of 108? to 10?5 Mo (Table 3), which corre- 
spond to a gas fraction log My,/M.. of 0.5-1096 with a median 
value of 2%. The two H 1 non-detections of FAST further con- 
strain the atomic-to-stellar mass ratio to 0.596, which is con- 
sistent with our expectations of 1%. The noise level of our 
new FAST observation is ~ 3—7 times better than that from 
the GBT or Arecibo, which is due to the large dish size thus 
high sensitivity of FAST compared to GBT; and significantly 
longer integration times relative to the ALFALFA survey. 
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Table 3. Upper Limits of H1 Mass 


Plate-IFU o log Mu tim Notes 
(mJy) (Mo) 
(1) (2) (3) (4) 
7977-3703! 2.6 9.0 2 
8081-6102 1.2 8.9 2 
8084-6103! 1.4 9.3 2 
0.2 8.2 1 
8728-3701? 1.3 8.9 2 
0.4 9.3 1 
8950-12705 2.6 8.9 2 
8952-12701 2.6 9.0 2 
9195-3702 0.3 8.8 1 
9512-3704 3.9 9.8 2 


Column (1): the MaNGA Plate-IFU. Columns (2): noise level of the 
Hı line profile. Column (3): upper limit of the Hı mass, assuming 
a 3c detection with a line width of 200 kms~!. Column (4): the H1 
spectrum is from 1: FAST and 2: H -MaNGA observation. 


' The Hı spectra of MaNGA 7977-3703 and MaNGA 8084-6103 
are classified as detections in H1-MaNGA. 

2 The Hr emission of MaNGA 8728-3701 is contaminated with 
that of its companion MaNGA 8728-12701, thus we calculate 
the upper limit of their total H1 mass. 


3.3 Hr Absorption in MaNGA 8155-6101 


Remarkably, MaNGA 8155-6101 has a ~ 6c: H1 absorption 
signal (Figure 4). MaNGA 8155-6101 is a merging galaxy 
system with a companion CGCG 391-011 at a velocity dif- 
ference of 360 km s^! and a projected physical separation of 
~ 30 kpc (44"). To investigate the absorption feature, we fol- 
low Wolfe & Burbidge [125] and Allison et al. [3], and the 
H1 column density from the absorption line signal is deter- 
mined as 


Ng, = 1.823 x 10 T. { 1(V)dV, 
18 T spin a 
~ 1.823 x 10 ES Tops (V)dV, 


where T(V) is the optical depth at a given velocity V, Tops(V) 
is the observed optical depth (panel e of Figure 4), Tspin is 
the mean harmonic spin temperature of the gas, and f is 
the covering factor of Hı absorption. The Hı spin temper- 
ature ranges from 10 to 1000 K by analysing the H remission 
and absorption profiles towards individual clouds located in 
front of Galactic and extra-galactic radio continuum sources 
[50, 77, 80, 108]. The covering factor f varies from 0.1 to 
1 [27] for galaxies with z <1, and we often assume it to be 
1 due to lack of information [78]. We adopt Tspin = 100 K 
and f = 1 following Carilli et al. [21]. This results in a Ht 
column density of 5.5 x 107° cm~?, which is within the typi- 
cal range of 102-10?! cm~? [80]. But given the uncertainty 
in the spin temperature and the covering fraction, the uncer- 
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tainty of Hı column density could be one or two orders of 
magnitude. The center of the absorption signal is precisely 
consistent with the optical central velocity of MaNGA 8155- 
6101, while the emission line signal on the blue-shifted side 
suggests a potential association with the neighboring com- 
panion, CGCG 391-011. 


Based on spatially resolved stellar kinematics, 
MaNGA 8155-6101 has a regular rotation pattern and high 
velocity dispersion within one effective radius (Panel c of 
Figure 4). However, the ionized gas kinematic traced by Ha 
emission shows significant perturbation. The CO emission 
line profile is quite extended with a line width ~ 600 km s^, 
even though the S/N is low [71]. The broad CO line width 
is rare considering there are only two galaxies with CO line 
widths larger than 600 km s^! among 333 detections in the 
xCOLD GASS sample [97], which may suggest that the CO 
in two interacting galaxies are perturbed and mixed. There- 
fore, the molecular gas in MaNGA 8155-6101 are potentially 
perturbed, either by AGN or galaxy interaction. Furthermore, 
its SFR, specific star formation rate (SSFR), and star forma- 
tion efficiency are typical for a massive quenched galaxy 
with the molecular gas fraction of 1% [71]. Thorp et al. 
[115] found that the central bursts of star formation in merg- 
ers require centralized enhancements in gas fraction, thus the 
broad CO line width and normal SFR in MaNGA 8155-6101 
may suggest that gas is distributed extensively thus failed to 
fuel central star burst activities. 


We employ three types of BPT diagrams ([N II], P1P2, 
and [S II]) to classify the ionization state based on the spa- 
tially resolved IFU data from MaNGA. MaNGA 8155-6101 
is classified as an AGN host based on BPT diagrams (panels 
i-n of Figure 4), and the fraction of AGN or Seyfert pixels de- 
creases from ~100% to 50% with the transition from [N II], 
P1P2, to [S II] BPT diagrams. The high ionization states is 
likely due to perturbations fuelling an AGN or driving shocks 
[60]. 


The possible overall scenario is that galaxy interaction dis- 
turbs the gas distribution in the galaxy outskirt and then re- 
sults in gas inflow and turbulent motion, which in turn trig- 
gers the AGN [46, 112, 121] within MaNGA 8155-6101. To 
gain further insights into the gas distribution and kinemat- 
ics under the effects of merger and AGN, future spatially- 
resolved H r observations may provide valuable details of the 
gas kinematics in the circum-galactic medium, intergalactic 
medium, and around the central black hole. 
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4 Atomic-to-Molecular Gas Conversion 


The atomic-to-molecular gas conversion is predominantly fa- 
cilitated in regions with high mid-plane pressure, weak ultra- 
violet radiation fields, and high metallicity, which holds both 
theoretically [41, 62] and observationally [66]. Intriguingly, 
a weak correlation between My, and My, in nearby galaxies 
suggests that the physical conditions governing the atomic- 
to-molecular gas conversion are non-trivial [23, 95]. To fur- 
ther constrain the conditions that influence the gas conver- 
sion, we carefully scrutinize whether the molecular-to-atomic 
gas ratio (Rmo = log My,/Mu,) depends on various physical 
properties. 


4.1 Stellar Mass Distribution and Star Formation 


We investigate the dependence of Rmo on M., Hs, NUV-r 
color, and sSFR in Figure 5 using the HICO-MaNGA sam- 
ple in comparison with the HICO-Spec sample (Figure 5). 
Galaxies with at least one detection in H1 or CO are shown. 
The NUV-r color is obtained from the NASA-Sloan Atlas 
(NSA) catalog [10], and its typical uncertainty is 0.2 dex 
[126]. The p value? and Kendall’s t coefficient indicate that 
Rmo is correlated with u., less tightly but still correlated with 
M., and NUV-r color, and not related to the sSFR. The best- 
fit linear relations in panels (a-c) of Figure 5 are derived from 
the orthogonal distance regression. The typical uncertainty of 
Hı and CO detection is 0.1 dex for log My, and 0.2 dex for 
log My, for our sample. While fitting the linear relations, we 
assume a typical uncertainty of 0.2 dex for log My, and 0.4 
dex for log My, for these non-detections. The uncertainties 
of p-value and Kendall's t coefficient are obtained by boot- 
strap realizations, which considers the uncertainties of both x 
and y axes. 

The HICO-Spec sample has a Rmo; ranges from -1.9 to 
0.8 with a median value of ~ —0.7 [23]. The values of Rag 
for our sample cover a similar range. The HICO-MaNGA 
sample shows a trend consistent with that of HICO-Spec and 
xCOLD GASS: Rq increases monotonically with the in- 
creasing of M., H+, and NUV-r color [15,16,23,94,95] with 
a scatter of ~0.5 dex. 

The positive relation between M. and Ryo! is consistent 
with literature studies. The typical value of logarithmic 
molecular-to-atomic gas ratio Rmo is 0.6 for S0/Sa galax- 
ies and —0.6 for Sd/Sm galaxies [93, 116, 129], which de- 
creases statistically from Sa to irregular galaxies (high mass 
to low mass, [15,130]). This weak relation holds true for both 
constant aco and luminosity-dependent aco [16]. Hydro- 
dynamical and semi-analytical simulations have also success- 
fully reproduced this relation [64, 87]. Therefore, the effi- 


3) Following convention, we consider two parameters are correlated if p « 0.05. 
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Figure 4 The physical properties of MaNGA 8155-6101, which is an interacting galaxy with H1 absorption signal. Panel (a) shows the SDSS g-band image 
of MaNGA 8155-6101 with its companion CGCG 391-011, and the field of view is 2.9’ x 2.9’ (same as FAST beam size). Panel (b) is the SDSS gri color 
composite image of MaNGA 8155-6101. Panel (c) and (d) are the velocity field and velocity dispersion of the stellar component for MaNGA 8155-6101. Panel 
(e) is e* as a function of velocity for MaNGA 8155-6101 and its companion, where r is the optical depth of the absorbing gas following Allison et al. [3]. 
Panels (f), (g), and (h) are the Ha image, velcoity field, and velocity dispersion field of MaNGA 8155-6101. Panels (i)-(n) represent the classification and 
spatial distribution of three ionization states distinguished by three different BPT diagrams. 


ciency of atomic-to-molecular gas conversion varies along 
different types of galaxies. 


We consider that the dependence of Rmoi on u+ is primar- 
ily due to the regulation of mid-plane pressure [41], which is 
within the expectation for a spatially-resolved view: Rinoi c 
uo (us + Us Vg /V)°S [62], assuming a thin disk of uniform gas 
and stars in balance [11]. Here, 4g is the gas surface density, 
and v, and v, are the vertical velocity dispersion of the gas 
and stellar components, respectively. In a spatially resolved 
view, the molecular-to-atomic gas surface density ratio is 
tightly proportional to u, observationally [66]. These results 
support the scenario that the atomic gas in high-density re- 
gions is more easily converted to molecular gas. This is also 
consistent with the observation that massive early-type galax- 
ies have higher values of Ryo! than late-type galaxies [15,94]. 


However, the HICO-MaNGA sample exhibits a steeper 
slope between Ry and u, compared to that of the HICO- 
Spec or xCOLD GASS sample. The reason may be that the 
HICO-MaNGA sample (1) mainly consists of star-forming 
galaxies, and (2) have log u, < 8.7 Mo kpc^?. In the bulge- 
dominated cases with high stellar mass surface densities u, > 
1057 Mo kpc^?, the cold gas reservoir can be efficiently con- 
verted to stars. Below this ju, threshold, higher detection rates 
are found in both molecular and atomic gas measurements 
[22, 95]. As galaxies become more bulge-dominated, both 
the atomic and molecular gas fractions decrease significantly 
as well as their detection rates [95]. The flattening of the rela- 
tions u, versus Rmo for the xCOLD GASS sample is evident 
above the characteristic thresholds of log pu. ~ 8.7 Mo kpc^?, 
which is mainly due to cold gas depletion in mostly quies- 
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Figure 5 The molecular-to-atomic gas ratio as a function of the log M, (panel a), log yz. (panel b), NUV—r color (panel c), and sSFR (panel d). In each 
panel, the black (detection) and grey (non-detection) symbols are galaxies from HICO-Spec, the red symbols are galaxies from the HICO-MaNGA sample. 
The Kendall’s 7 coefficient and p value are shown in the upper left corner (grey: galaxies in HICO-Spec with at least one detection in H 1 or CO; black: galaxies 
in HICO-Spec with both detections in Hı and CO; red: HICO-MaNGA). The Hı or CO non-detection is shown as up and down arrows. The best-fit linear 
relations based on the orthogonal distance regression are shown as solid lines in panels (a)-(c) with the shaded region denoting the 1c scatter of the fit. 


cent, early-type galaxies [94]. The HICO-MaNGA sample 
is mainly star-forming galaxies and not depleted in cold gas, 
thus the flattened trend is not observed. 


We investigate the relation between Ry; and star forma- 
tion activity in the bottom panels of Figure 5, which show a 
weak and no relation, respectively. NUV-r color is an in- 
dicator of star formation history or sSFR in galaxies [126]. 
Saintonge et al. [95] also showed Ry is an increasing func- 
tion of NUV-r, but this relation may be a direct result of 
the strong relation between atomic gas fraction and NUV-r 
[22, 23, 136]. The NUV-r color is sensitive to star formation 
within several 100 Myrs [65], but sSFR is more sensitive to 
star formation for a longer time scale (several Gyrs), because 
it quantifies the current star formation activity relative to the 
existing stellar mass. So the tighter relation between Rmol 
and NUV-r color compared with that of sSFR may suggest 
that the gas conversion is related with the star formation in 
relatively short timescales. 


The HICO-MaNGA and xCOLD GASS samples show no 
relation between Ryo and sSFR (panel d of Figure 5). The 


weak relation is consistent with the results based on 18 detec- 
tion out of 33 galaxies from the ALMaQUEST sample based 
on Hr-MaNGA data release 1 [71]. From late-type to spi- 
ral galaxies (SSFR> 107!? yr7'), the value of Ry increases 
slowly and then becomes constant in spirals [15]. The global 
SSFR is related to the atomic gas fraction (e.g., [35,51, 138]) 
and molecular gas fraction [15] with large scatters. For galax- 
ies on the star-forming main sequence, their variations of Ryo 
are mainly influenced by their H 1 reservoirs [23], but the star 
formation activity relies more on H, content [57]. In addi- 
tion to the total gas content, the central concentration of cold 
gas within the optical disk may be more directly related to 
star formation [122, 132], especially when there is an abun- 
dant gas reservoir. Our results suggest that the global conver- 
sion from H1 to H» does not significantly influence the global 
sSFR in non-starburst galaxies. 


4.20 Cold Gas Distribution 
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Figure 6 The molecular-to-atomic gas ratio as a function of the H1 profile asymmetry Ar (panel a) and profile shape K (panel b) for the HICO-Spec sample 
with H 1 detection. The black dots are CO detection, and the grey arrows are CO non-detection. The blue circles are the median trend with the standard deviation 
of median shown as the errorbar, and the blue shaded regions are the 25th and 75th percentiles of the distribution. The p value and Kendall v are shown in the 


upper left corner in each panel. 


Even though the spatially resolved H 1 and CO are not simul- 
taneously available for galaxies with a large sample size (e.g., 
[66, 135]), the asymmetry and relative concentration of cold 
gas could be inferred from their global profiles [131, 132]. 
Meanwhile, spatially resolved gas asymmetry is reflected in 
asymmetry parameters derived from the global profile [89]. 
The double-horned profile with a deep central trough may 
suggest centrally depressed H1 distribution [92]. So the 
global profile encodes the gas spatial distribution. 

We investigate the dependence of Ray on the Hı pro- 
file asymmetry Argi and profile shape Kg; for galaxies with 
H1 detection in the xCOLD GASS and HICO-Spec samples 
(Figure 6). Our analyses indicate a weak or no relation be- 
tween Ry and the atomic gas distribution, inferred from ei- 
ther asymmetry Ap pı or profile shape Kyr. The p value for 
Rmoi Versus Ap yy is slightly larger than 0.05, thus their re- 
lation is not significant. The p value for Rmo: versus Kyy is 
significantly smaller than 0.05, which suggest that the gas 
conversion may be related to the H1 profile shape K. 

Asymmetric features in the distribution and kinematics of 
Hı can be attributed to either internal perturbations, such as 
stellar [26] and active galactic nuclei feedback [79], or ex- 
ternal perturbations, such as accretion [102], stripping [123], 
and merging [13]. We caution that the profile asymmetry of 
massive mergers is not higher than that of the control galax- 
ies after correcting for the effects of S/N [139]. Based on 
high-resolution, spatially-resolved H 1 distributions in nearby 
galaxies, Reynolds et al. [89] argued that gas-removal mech- 
anisms, such as tidal interactions and ram-pressure stripping, 
are the most probable mechanisms for generating asymmet- 
ric H 1 distributions. The weak relation between Ryo} and H t 
asymmetry may suggest that the gas-removal mechanisms do 


not significantly influence the H r-to-H» conversion. 

The panel (b) of Figure 6 demonstrates that galaxies with 
higher Rmo tend to exhibit more single-peaked H i profiles in 
our study, suggesting that H ris relatively more centrally con- 
centrated within the optical disk. Perturbations may drive H1 
inflow [34, 37, 105] and promote the atomic to molecular gas 
conversion [110]. However, we note that the inclination an- 
gle, stellar mass, optical concentration, and gas mass cannot 
be well controlled in this work as done in Yu, Ho, & Wang 
[132] to better infer relative gas concentration within the op- 
tical disk for nearby galaxies due to the limited sample size. 

Hı profile shape asymmetry Ac pı shows similar results 
as that of Ap gr, and profile concentration Cy py also demon- 
strates comparable trends to Ky. Therefore, the H1 distribu- 
tion does not significantly regulate the atomic-to-molecular 
gas conversion. 

Figure 7 illustrates the trends of Ry as a function of CO 
profile asymmetry Apo and profile shape Kco. The uncer- 
tainty of CO profile asymmetry is high, because there are lots 
of marginal detections with low S/N. Same as Figure 6, the 
dependence of Ryo; on CO distribution is weak or no with 
p=0.131 for Rmo; versus Agco and p=0.031 for Rmo versus 
Kco. The CO redistribution or asymmetry can be caused by 
the AGN activity [25, 113] or bar structures [109], so the gas 
conversion does not influence the molecular gas distribution. 

Even though the trend is as weak as that of H 1, it is note- 
worthy that the trends in CO are opposite to those in H1. The 
contrasting trends suggest that the distribution and kinematics 
of Hı and CO in galaxies are not identical, which is consis- 
tent with previous findings by comparing global profiles of 
Hı and CO [30,90]. The spatially resolved H r and CO asym- 
metry is weakly correlated for galaxies in the Virgo cluster, 
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Figure 7 The molecular-to-atomic gas ratio as a function of the CO profile asymmetry Ar (panel a) and profile shape K (panel b) for the xCOLD GASS 
sample with CO detection. The black dots are H 1 detection, the grey arrows are H 1 non-detection, and the blue circles shows the median trend. The standard 
deviation of median is represented by its errorbar, and the blue shaded regions are the 25th and 75th percentiles of the distribution. The p value and Kendall t 


are shown in the upper left corner in each panel. 


but a stronger correlation for the galaxies strongly perturbed 
by environmental effects [90]. Specifically, H1 disks tend 
to be more fragile and extended (such as [31]), but CO is 
more centrally concentrated within the optical disk. Using 
global profiles in the xCOLD GASS survey, we find no cor- 
relation between the global profile asymmetry and shape of 
Hı and CO, which indicates that the distribution of atomic 
and molecular gas in nearby galaxies is not tightly linked. A 
comprehensive analysis of global H 1 and CO profiles will be 
fully investigated in an upcoming paper. 


4.3 Ionization States 


We investigate the dependence of the atomic gas fraction log 
(My,/M,.) and molecular gas fraction log (My,/M.) on the 
fraction of different ionization states for the HICO-MaNGA 
sample employing the P1P2 BPT diagram (panel a and b of 
Figure 8). We have found that the H 1 region fractions mono- 
tonically declines to ~ 50% with decreasing atomic gas frac- 
tion. A similar trend is observed for the molecular gas frac- 
tion. As the fraction of atomic or molecular gas increases, the 
proportion of H n regions in galaxies increases from ~40% to 
100%. Galaxies with log (My,/M,)= —0.8 tend to have Hu 
region fractions of greater than or equal to 80% (panel a of 
Figure 8). Thus the most gas-rich galaxies are predominantly 
star-forming. A similar trend appears in the molecular gas 
fraction. These trends are consistent across the three types 
of BPT diagrams: nearly identical results for the [N II] and 
P1P2 BPT diagrams, whereas the [S II] BPT diagram exhibits 
a lower fraction of composite regions (<20%) and AGN re- 
gions (X596). 

To disentangle the dependence of log (My,/M..) and log 


(My,/M.) on M., we investigate the relations of the verti- 
cal offsets with respect to the median trends as a function 
of M, and the fractions of three components (panel c and 
d of Figure 8). With the increasing of A log (My,/M,) and 
A log (My,/M.,), the fraction of emissions from the Hm re- 
gions increases monotonically, similar to the trends shown in 
panels (a) and (b). Therefore, the positively-correlated trend 
indicates that an abundant cold gas reservoir, both atomic 
and molecular, promotes wide-spread star formation within 
galaxies. 

We investigate how the molecular-to-atomic gas ratio 
(Rao) and the galaxy ionization states influence each other 
(Figure 9). All galaxies of HICO-MaNGA have Ha equiv- 
alent widths larger than 1 A, both within the field of view 
and within one effective radius, thus none of them is emis- 
sion line-less galaxy [7]. We adopt the mask of each emis- 
sion line and require the S/N of each spaxel higher than 10 
to minimize the effects of noise [124]. There are 3 galaxies 
(MaNGA 8655-12705, MaNGA 12700-12702, and MaNGA 
12769-6104) have numbers of spaxel smaller than 100 after 
adopting the S/N threshold, which may lead to high uncer- 
tainties of the fraction of different components, thus we re- 
move them in the following discussion. 

With the increasing of Rmo from —1 to 0, the H u region 
fractions slightly decreases from 10096 to 8096 (panel a of 
Figure 9). The decreasing of the fraction of the Hm regions 
may be caused by the decreasing of cold gas fraction (see Fig- 
ure 8). Because the decreasing trend has a large scatter and 
the BPT classification is a mixed effect of metallicity and ion- 
ization states, we divide the HICO-MaNGA sample into three 
Rmo bins with same galaxy numbers and investigate their dis- 
tribution in the P1P2 BPT diagram (panels b-d in Figure 9). 
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Figure8 The fraction of different optical emission line components classified by PIP2 BPT diagram as a function of the (a) atomic gas fraction, (b) molecular 
gas fraction, (c) deviation of the atomic gas fraction, and (d) deviation of the molecular gas fraction. In each panel, the small dots show the fraction of each 
individual galaxy, the squares show the median value of the x and y-axis, and the shaded regions represent the 25th and 75th percentiles of the distribution. 


The galaxy ionization state increases with the increasing of 
P1 parameter, and the metallicity decreases with the increas- 
ing of P2 parameter [52]. With the increasing of Rmo, P2 
decreases significantly and P1 evolves into two branches. 

From Rmo < -0.75 to Rma > —O0.75, the fraction of 
composite regions increases by 23-696. Efficient atomic-to- 
molecular gas conversion will promote star formation and en- 
hance stellar feedback. The increasing of composite region 
may be the mixed effects of star formation, shock excitation, 
and/or AGN activity [60]. On the other hand, the values of 
P1 decrease significantly with the increasing of Ryo, but the 
regions are still classified as H n-dominated emission. The P2 
values are anti-correlated with galaxy metallicity, thus higher 
values of Rmo corresponds to higher galaxy metallicity. More 
researches are vital to reveal detailed physical mechanisms 
behind the bilateral evolution of P1. 

The efficiency of atomic-to-molecular gas conversion has 
been found decreasing monotonically with with the increas- 
ing of galactic radius [66], thus we evaluate the dependence 
of global Rmo: on the fraction of different ionization states 
within a given galactic radius range in our samples. We 
divide spatially resolved BPT diagrams (see the example 
in Figure 4) into two or three radius ranges: R € O.5R,, 
O05R, < R <1.5R., L.58R, < R x 2.R,, where R, is 
the effective radius of each galaxy. The results in Figure 10 
show that with the increasing of galaxy radius, the fraction of 


non-star-forming regions decreases from ~ 20% to 2%. It is 
worth noting that more than 80% of our sample have an AGN 
fraction below 5% in their BPT diagrams, and more than half 
of our sample have a fraction of the combination of AGN and 
composite regions lower than 10%, thus our sample exhibits 
a low fraction of non-star-forming regions. The non-star- 
forming emission could be mainly caused by AGN or shock 
activity [60]. The HICO-MaNGA sample contains 15 galax- 
ies with a weighted AGN and composite fraction higher than 
15%, which tend to be AGN hosts [127]. However, only one 
galaxy, MaNGA 7977-9101, has a typical Ha velocity disper- 
sion higher than 150 km/s, which could be caused by galactic 
wind shocks or AGN feedback [29, 60]. Thus, it is difficult to 
distinguish the main physical drivers for the decrease in the 
fraction of star-forming regions with the increasing of Rao 
using the current data. 

We specifically examined 5 galaxies characterized by a 
high molecular-to-atomic gas ratio and a low fraction of Hu 
regions (Rmo > — 0.3 and the fraction of Hu regions < 
50%): GASS 4030, GASS 11071, MaNGA 8081-12703, 
MaNGA 9194-3702, and MaNGA 8655-3701. These galax- 
ies exhibit distinct characteristics: bright Ha core and an old 
stellar population (high values of D,,4000) surrounding the 
Ha core. Their main disks show significant AGN or com- 
posite emission. These can be shown in Figure 11, where 
we show the images of the optical, D,4000, Ha, and BPT of 
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Figure 9 Panel (a): The fraction of different optical emission line components classified by P1P2 BPT diagram as a function of the molecular-to-atomic gas 
ratio (Rmoi). The classification and spatial distribution of three ionization states distinguished by the P1P2 BPT diagram for galaxies with (b) Roi € —0.75, (c) 
—0.75 < Rmol < —0.40, and (d) < Ry = —0.40. In panels (b)-(d), the gray scale and contours show the density distribution for all spaxels in each panel, and 
the red cross shows the density peak of the distribution in panel (b) to highlight the difference of the distribution in three panels. The galaxy numbers in each 
panels (b)- (d) are the same, and the fractions of Hm composite, and AGN-dominated regions are shown. 
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Figure 10 The fraction of different components classified by BPT diagram as a function of the molecular-to-atomic gas ratio. From left to right, the galaxy 
radius range is (a) R € 0.5R,, (D) 0.5Re < R € L5R,, and (c) 1.5Re < R € 2.R,. In each panel, the small dots show the fraction of each individual galaxy, 
the squares show the median value of the x-axis and y-axis, and the shaded regions are the 25th and 75th percentiles of the distribution. 


MaNGA 8655-3701: a typical galaxy among these outliers. 
Both the stellar and Ha components display a rotation pat- 
tern, and its bright central core has high velocity dispersion. 
The relative distribution of H 1 and H5 may be offset, because 
there is significant difference between their global profiles. 


Consequently, we found that a high global value of Ry 
is not necessarily linked to a high fraction of Hm regions. 
They are separate processes within the baryon cycle in certain 
cases. The decline in Hu region fraction at high molecular- 


to-atomic ratios can be mainly attributed to high velocity dis- 
persion (see the Ha velocity field in Figure 11) in the inter- 
mediate region (0.5R, < R < 1.5R,) and extended cold gas 
distribution. These factors hinder the connection between H r, 
Hp, and star formation. A similar effect is found that inten- 
sive star formation may promote gas velocity dispersion, but 
nonaxisymmetrical torques can prevent the gas from being 
gravitationally unstable [61]. 
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Figure 11 The physical properties of MaNGA 8655-3701: panel (a) shows the SDSS composite image, panel (b) and (c) are the velocity field and velocity 
dispersion of the stellar component, panel (d) is image of D,4000, panels (e)-(g) show the Ha image, velocity field, and velocity dispersion, and panel (h) 
presents the spatially resolved BPT diagram. In panel (h), the red, blue, and green spaxels showing the AGN, composite, and H i regions, respectively. 
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Figure 12 The molecular-to-atomic gas ratio as a function of the central oxygen abundance (panel a), the oxygen abundance at one effective radius (panel b), 
and the slope of oxygen abundance within 0.5 to 2 effective radius (panel c). The oxygen abundance and its slope are derived from O3N2 calibrator taken from 
Pipe3D [101]. In each panel, the p value and Kendall coefficient t as well as their uncertainties are shown in the upper left corner. In panels (a) and (b), the 
best-fit linear relations based on the Orthogonal distance regression are shown as the black line, with the black shaded region denoting the 10 scatter of the fit. 


4.4 Metallicity 


Metallicity acts as a catalyst in the process of converting Hi 
Gas 
metallicity serves as a direct indicator of dust [33], and ef- 


to Hz gas, which enhances gas conversion [41, 62]. 


ficient Hz formation occurs in the surface of dust grains [45]. 
So metallicity plays a positive role in gas conversion. Boselli 
et al. [15] collected optical emission lines and utilized the 
O3N2 calibration following Pettini & Pagel [84]. Their re- 
sults indicated a correlation between Ryo and the oxygen 


abundance 12+log (O/H) globally. However, the concern 


about the non-uniformly measured oxygen abundance due to 
the limitation of the quality of the data could affect the con- 
clusions. 


To address this, we re-investigate the dependence of Rmol 
on the oxygen abundance 12 + log(O/H) by using the data 
from the MaNGA IFU. The oxygen abundance and its slope 
are taken from Sánchez et al. [101], which is determined 
using the O3N2 = ([O IIIJ/H8)/[N II]/Ho) calibration dis- 
cussed above [84] for star-forming regions classified by the 
[N II] BPT diagrams with with EW(Ha) > 3À. The slope of 
the oxygen abundance is derived from the fitting results ob- 
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tained between 0.5 and 2.0 Re. As shown in Figure 12, the 
conversion from atomic to molecular gas is enhanced with 
the increasing of the oxygen abundance at central and one 
effective radius. The oxygen abundance calibrated by the 
O3N2 index is widely used in high-metallicity regime at so- 
lar and super-solar metallicity regime [74, 84], where [N II] 
saturates. The oxygen abundance based on the R23 = ([O 
II]--[O III]J/H8) and N2 = [N II]/Ha calibration [101] returns 
to similar results. 

So metallicity plays a positive role in gas conversion. 
Compared to the central metallicity, the metallicity at one 
effective radius shows a stronger relation with global Ryo. 
It suggests that gas conversion predominantly occurs within 
the optical disk, even though the central region exhibits 
higher conversion efficiency. We also found that the xCOLD 
GASS sample represents a positive relation between Rio} 
and global gas-phase metallicity. 

Our data show a weak or no positive relation between Ry 
and the slope of the oxygen abundance: p=0.112 and t=0.11. 
The positive slopes of oxygen abundance is likely the result of 
gas accretion [118], which dilutes metallicity. However, the 
gradient of oxygen abundance is slightly dependent on stel- 
lar mass [8,119,120], although EAGLE simulations show no 
clear trend between them [117]. The physical drivers of the 
observed metallicity gradients are complicated, including the 
radial variations in the SFR [85, 106], gas motions [82, 88], 
and IMF [47, 75]. Therefore, the lack of correlation between 
Rmo and the oxygen abundance gradient does not necessarily 
rule out the effects of gas accretion. 

The Kendall’s r coefficient for global Ry versus the oxy- 
gen abundance at one effective radius is higher than other re- 
lations in Figure 5, 6, 7, and 12. Among the dependencies 
studied in this work, metallicity within the optical disk may 
be the most important factor regulating the H 1-to-H, conver- 
sion in galaxies, which is consistent with literature studies 
[36, 45]. 


5 Summary 


The atomic-to-molecular gas conversion is an important step 
of the galaxy baryon cycle. We obtained the H1 spectra line 
observations of FAST for the ALMaQUEST sample. We 
derived 42 emission lines, 1 absorption lines, and 8 non- 
detection after combining data from Arecibo and GBT. Our 
FAST observations are with a noise level of approximately 
0.2 mJy with a velocity resolution of ~ 6 km s^!, which is 
much better than that of H -MaNGA. Thanks to the high sen- 
sitivity of FAST, the H1 detection rate is 80%, rendering it an 
ideal instrument for extragalactic H I researches. 
Additionally, we compile H 1 and CO spectra or data from 
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the xCOLD GASS, ALLSMOG, and JINGLE surveys to con- 
struct two samples: HICO-Spec and HICO-MaNGA, which 
contains around 300 and 100 galaxies, respectively. The 
HICO-MaNGA sample has H t, CO, and optical IFU obser- 
vations. This sample is predominantly composed of star- 
forming galaxies with a stellar mass range of 10??—10!!? Mo 
and a redshift range of 0.02 x z x 0.06. We analyzed the 
H1 and CO spectra using our recently developed “curve of 
growth" method [131, 133] to measure parameters such as 
the total flux F, profile asymmetry Ar, and profile shape K. 
The main results are briefly summarized as below. 


e The molecular-to-atomic gas ratio Ryo) is positively re- 
lated to the M., u., and NUV-r color, but shows no dis- 
cernible relation with the sSFR. In contrast to the results of 
Saintonge & Catinella [94], our sample does not show a flat- 
tened trend in the relation of u, and Rmo. This difference 
may be attributed to the star-forming nature of our sample 
and low stellar mass surface density (log yu. < 8.7 Mo kpc^?). 
Gas conversion depends on mid-plane pressure, which is pri- 
marily proportional to surface density, assuming a thin disk 
with uniform gas and star distribution [11]. 


e The relationship between the Rmo and the distribution of 
H1and CO is weak, while the gas distribution is inferred from 
the asymmetry and shape of their global profiles. Galaxies 
exhibiting higher Hı asymmetry or more single-peaked Hı 
profiles tend to have higher values of Riot, while the opposite 
trend is observed for CO. Because the relation is weak, we 
cannot draw definitive conclusions regarding the relationship 
between the cold gas distribution and H r-to-H» conversion. 


e With the increasing of Ry, the fraction of H u region 
decreases within 1.5R,, likely due to the influence of shocks. 
We determine the spatially resolved ionization states using 
the BPT diagram introduced by Ji & Yan [52] and examine 
the dependence of Ry) on the fractions of different regions. 
The fraction of Hm regions increases with the increasing of 
atomic and molecular gas fraction, which is still evident after 
controlling the effects of stellar mass. Consequently, cold gas 
reservoirs enhance star formation within galaxies. 


e However, the fraction of Hm regions decreases with the 
increasing of Rag. The decreasing of H u regions fraction is 
related to the decreasing of cold gas fraction. With the in- 
creasing of Rmo1, the metallicity increases and the ionization 
states diverge. The fraction of composite regions increases by 
2%, which may be due to the stellar feedback. Therefore, a 
high molecular-to-atomic gas ratio does not necessarily lead 
to a high fraction of H u regions, as this is also regulated by 
galaxy kinematics and gas distribution. 


e Furthermore, galaxies with higher oxygen abundance 
12 + log (O/H) tend to have higher values of Ryo, indi- 
cating that metellicity acts as a catalyst enhancing atomic- 
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to-molecular gas conversion. The oxygen abundance at R, 
proves to be more critical than the central oxygen abundance 
or oxygen abundance gradients. The slope of oxygen abun- 
dance is not tightly related with the global Ryo, which may 
suggest that gas accretion does not significantly promote the 
atomic-to-molecular gas conversion. 


In summary, the atomic-to-molecular gas conversion 
mainly depends on stellar mass surface density and metal- 
licity. The most efficient atomic-to-molecular gas conversion 
occurs in massive or metallicity-rich galaxies. Over a longer 
timescale, the feedback mechanisms may play a role in the 
galaxy ionization states, thereby influencing the gas conver- 
sion. To gain deeper insights, multi-wavelength observations 
for stars, ionized gas, atomic gas, molecular gas, and dust 
with similar resolutions are needed to further constrain the 
atomic-to-molecular gas conversion in nearby galaxies. 
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